Nonlinear a model approach for level correlations in chiral disordered systems 
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We study level correlations of disordered systems with chiral unitary symmetry ( AIII symmetry) . 
We use a random matrix model with a finite correlation length to derive a supersymmetric nonlinear 
0" model. The result is compared with existing results based on other models. Using the methods 
by Kravtsov and Mirlin (Pis'ma Zh. Eksp. Teor. Fiz. 60, 645 (1994) [JETP Lett. 60, 656 (1994)]) 
and Andreev and Altshuler [Phys. Rev. Lett. 75, 902 (1995)], we calculate the density of states 
and two-level correlation function. The result is expressed using the spectral determinant as in 
traditional nonchiral systems. We discuss the renormalization of the mean level spacing which is 
not present in the traditional systems. 

PACS numbers: 05.30.-d, 05.45.-a, 72.15.Rn, 71.30.+h 
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I. INTRODUCTION 

The classification of disordered systems is based on 
symmetries of the Hamiltonian. According to invariance 
properties under time-reversal and spin rotation, three 
symmetry classes - unitary, orthogonal, and symplectic 
- are well known since the work by Wigner and Dyson 
The modern classification is based on the notion 
of symmetric spaces and indicates that ten univer- 
sality classes exist. Although there was an early effort 
at a universality classification in 80's 3], the additional 
seven classes did not attract much attention until phys- 
ical applications were found 0, |^. The importance of 
chiral symmetry in disordered systems was first noticed 
in Ref.j^] by using random matrix theory (RMT) in the 
context of quantum chromo dynamics (QCD) and meso- 
scopic quantum wires. In systems with chiral symmetry, 
eigenvalues of the Hamiltonian appear in pairs ±e and 
the origin e = plays a special role for level correlations. 

In order to analyze such systems, the supersymmetry 
method is known to be a useful tool for both pertur- 
bative and nonperturbative calculations. This method 
allows one to obtain a nonlinear a model with super- 
matrix fields as effective modes. One can discuss weak 
localization effects using perturbation theory, where an 
expansion in terms of diffusion propagators H29|) is per- 
formed. A diagrammatical interpretation is thus possi- 
ble, and weak localization implies a large conductance 
3^1, where g is proportional to the diffusion constant 
in the propagator. The localization property can also be 
discussed using the renormalization group method. This 
expansion is justified only for nonzero modes q ^ Q in 
the propagator. The zero mode sector contains a totally 
different contribution and gives the ergodic result g = oo. 
Using the zero mode, we can calculate level correlation 
functions scaled in terms of the mean level spacing Q. 
The result is nonperturbative, parameter-free, and uni- 
versal. We know that treating the zero mode perturba- 
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lively gives only the asymptotic form of the exact result. 

Thus it is important to notice the different roles of the 
zero and nonzero modes. At finite g, the nonzero modes 
modify the universal result of level correlation functions 
0, 0, • Kravtsov and Mirlin (KM) treated the zero and 
nonzero modes separately and found finite-g corrections 
to the universal result Due to technical problems, 
the result was restricted to the domain z <^ g where z is 
the scaled energy variable. Using another method, An- 
dreev and Altshuler (AA) considered the domain z ^ 1 
where the perturbative expansion makes sense 0. They 
reached the nonperturbative regime by noticing the ex- 
istence of a set of nontrivial saddle points. Considering 
the expansion around two saddle points the result was 
expressed using the determinant of the diffusion prop- 
agator, which is called the spectral determinant in the 
literature. Although their method did not treat the zero 
and nonzero modes separately, it was shown in Ref.|lCl|| 
that the separation, just as in KM's method, gives the 
same result. 

Using the derived result, the authors in Ref.|3 found a 
smearing of the singularity at the Heisenberg time in the 
form factor (the Fourier transform of the two-level cor- 
relation function). Furthermore, the use of the spectral 
determinant represents a link from disordered to chaotic 
systems. The authors in Rcf. 11] noticed that a similar 
treatment can be applied to general chaotic systems just 
by replacing the diffusion operator in the spectral deter- 
minant by the Perron-Frobenius operator. For a chaotic 
system, the expression of the determinant using the trace 
formula was discussed in Ref . [l^ . Thus the expression 
using the spectral determinant is important for a uni- 
fied treatment of disordered and chaotic systems. The 
result was applied to critical statistics and the rela- 
tion to the density-density correlation in the Calogero- 
Sutherland model at finite temperature was discussed. 

In this paper we consider systems with chiral unitary 
symmetry. Starting from a chiral random matrix model 
with a finite correlation length, we derive a nonlinear 
a model and calculate the density of states (DOS) and 
two- level correlation function (TLCF). Our aim in this 
paper is not to discuss a specific model but to discuss the 
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generic properties of chiral symmetric systems. Actually 
the a model we use in this paper is believed to be appli- 
cable to a broad range of physical systems and we discuss 
the relation to other a models for specific systems. Then, 
we calculate the DOS and TLCF using a nonperturba- 
tive method which is equivalent to both methods by KM 
and AA. Out method is similar to that in Ref.To] and 
the zero and nonzero modes are separated explicitly. For 
chiral symmetric systems, the calculation using the KM 
method has been carried out in Ref . . In contrast to 
the approach in Refs.|3. IT^. we integrate the zero mode 
first, and then, treat the nonzero modes perturbatively. 
The advantage of this method is that all domains are 
treated in a unified way. We also discuss the effect of the 
DOS renormalization, which is specific for nonstandard 
symmetry systems. We restrict our discussion to chiral 
unitary symmetry (AIII symmetry) and the extentions to 
other chiral symmetric classes, chiral orthogonal (BDI) 
and chiral symplectic (CII), will be discussed elsewhere. 

The organization of this paper is as follows. In SeclTll 
starting from the random Hamiltonian, we derive the su- 
persymmetric nonlinear a model. It differs from the tra- 
ditional a model written in terms of a supermatrix Q by 
symmetries of the matrix and the presence of an addi- 
tional term. We discuss relations to other models. Next, 
the DOS and TLCF are calculated in SecHTIl In SeclTTI 
we discuss the effect of the additional term and the DOS 
renormalization. Sec^is devoted to discussions and con- 
clusions. 



II. SUPERSYMMETRIC NONLINEAR a 
MODEL 

A. Derivation 

In this paper we treat the Hamiltonian in the form 



H 



W 




(1) 



where W is an arbitrary matrix. This Hamiltonian pos- 
sesses chiral symmetry, which means that the eigenvalues 
appear in pairs ±6^. The matrix W can be a rectangu- 
lar matrix (n x to) as well as a square one (n — to). 
V — \n — m\ is the topological number and is equal to 
the number of zero eigenvalues of the Hamiltonian. Here 
we consider 1^ = and the extention to a finite v will be 
discussed elsewhere. Making a unitary transformation, 
we have 



H = 



-rii 



(2) 



fii 2 are n x n Hermitian matrices. Treating these ma- 
trices as random ones, we can obtain the original chiral 
RMT 4J. Due to the chiral structure of the Hamiltonian, 
two random matrices couple in the single Hamiltonian 
and nontrivial correlations of the single Green function 



are expected. We restrict our discussion to the chiral uni- 
tary ensemble, which means f2i_2 are arbitrary Hermitian 
matrices. 

We consider a system written in field theoretical form 

as 



I xy 



(3) 



where is the fermionic field operator and — J df^x. 
The random Hamiltonian H{x,y) has the chiral structure 
(0) and 



^lA^iV) = uji,2ix,y)ai\x - y\). 



(4) 



uji,2 are random matrices and are averaged using the 
Gaussian integral 

(•••) = [ Vuji^2i- ■ ■)exp [ {LJi{x,y)uji{y,x) 

J L ^ Jxv 



+uj2{x,y)uj2iy,x)) 



(5) 



where A is a free parameter. The function a{r) represents 
a finite correlation of the Hamiltonian. We assume the 
range of the correlation, denoted by rp, is large so that 
the saddle point approximation is applicable in the fol- 
lowing calculation. In the limit r ^ tq, a(r) ^ 1 and we 
have the fully Gaussian correlation. In the opposite limit 
r 3> To, we assume the correlation decays fast enough, 
e.g., a{r) ~ exp(-|r|/ro). 

This finite-range model is more realistic than chiral 
RMT in which all the matrix elements correlate with 
each other in the same way. The finite-range effect can 
be realized as a weak localization correction and a new 
energy scale Ec = D/L^ (Thouless energy), where D is 
the diffusion constant and L the system length, comes 
into the analysis. Another interesting situation is when 
the decay of the matrix is power-law. For a certain range 
of parameters this model reproduces the physics of the 
Anderson transition 15;]. Extentions of the present work 
to the power-law case are discussed in Ref.llq| . 

We mention related works [TtI IT^ ITgl I20ll2lj in which 
similar nonlinear a models were considered for systems 
with chiral symmetry. Our model is a simple general- 
ization of models used in 0, In other works, the 
random flux model ^^li the random gauge field model 
|2l|. and the partially quenched chiral perturbation the- 
ory as the low-energy model of QCD 19] were consid- 
ered. The derived nonlinear a models differ from the 
standard diffusion model for nonchiral systems by sym- 
metries of the matrix. Furthermore, an additional term 
was found in Refs. [T^ ITM Il9| although it was not found in 
Refs.pol l2l| . Here we rederive the a model and discuss 
relations to these models. In fact the additional term 
can exist and can be derived by a careful treatment of 
the massive modes integration. Although these models 
are different, we expect common low-energy properties. 
Our goal is to investigate them in the framework of the 
nonlinear a model. 
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Let us derive the nonlinear <j model using the super- 
symmetry method. Our derivation is similar to that in 
Refs.j3, 113' We first define the generating function 
for the single Green function. Following the Efetov's 
notation and conventions j6||, we define it as Zi[J] = 
/ 'D{^lJ, ip) exp(— £), with 



C 



-I I il}{x)[e^ 5{x — y) 

I xy 

-H{x,y) + kJ{x)5{x^y)]^{y), (6) 



where k = diag(l, —1) operates in superspace, -0 is a four- 
component supervector and •ip = ip'f . The source field J 
is a 2 X 2 matrix in chiral space. We take the ensemble 
averaging to obtain 



A{x,y)sti- p{x)p{y) 



xy 



-i / Tp{x) [e+ + kJ{x)] il){x), 

J X 

where A{x, y) = a^{x — y), and 

[p{x) - 'Ej;p{x)T,x] , 



(7) 



p{x) 



1 

72 



p{x) = syV(a;)V^(a;)Si/2. 



(8) 



Sx.z are the Pauli matrices in chiral space. The Hubbard- 
Stratonovich field Q is introduced in the standard way. 
After integrations over ^p and ip, we have {Zi[J]) = 
J VQexp{-Fi[J]) with 



Fi[J] 



2 



iA-')ix,y)stTQix)Q{y) 



str In e+E^ + JkY,^ + iX\/ AqQ 



(9) 



where Aq — A{x, y) ~ Tq. Q is a 4x4 supermatrix and 
has the same symmetry as p{x), which gives the condition 
{O,EJ = 0. ' 

We consider the saddle-point approximation. We are 
interested in the vicinity of the origin e = where chiral 
symmetry becomes important. At this point, the saddle 
point equation gives = 1 and the saddle-point mani- 
fold is obtained a.sQ = TTi^T where T is the inverse of T. 
Symmetries of the T-matrix were considered in Ref. |23|| 
and the explicit parametrization was obtained as 



T = y/l - P2 -ip_ p 




















ib ) 



(10) 



where a, h are real variables and cr, p grassmann vari- 
ables. In addition, we must take into account the mas- 
sive degrees of freedom which are not on the saddle point 
manifold. Usually, in nonchiral systems, integrations of 
the massive degrees of freedom do not give any contribu- 
tion. However, in the present case, the integrations give 



additional contributions written in terms of the massless 
modes. We can write the Q-matrix as Q = T{Yjz + 5Q)T 
where 5Q denotes the massive modes and changes the 
saddle point. Since the Q-matrix anticommutes with E^^, 
the structure oi 5Q in chiral space is determined as 



5Q 



Sq 
-5q 



(11) 



where 5g is a 2x2 supermatrix. This Q is substituted to 
the generating function and the functional Fi is expanded 
in powers of SQ. We have 



Fi[J] ^ f[°\j] + F^"" + F, 



(0) 



(12) 



where 



(0) 



- ' 2 



xy 

lire 

'WK 

it: 



F, 



(0) 



i?(a;,y)strQ(a;)Q(y) 
strE^Q(a;) 
str J{x)kTizQ{x), 
[{A-^){x.,y)+5{x-y)A-^'] 



(I) 



R{x,y) 



xy 



2\/A 
Ao 
2 

xstr SQ{x)5Q{y), 
Ao 
2 

xstr {2 [Tiy)Qix)Tiy) - E,] 6Qiy) 
+T{x)5Q{x)f{x)T{y)5Q{y)f{y) 
-5Q{x)5Q{y) + ---). (13) 

Q{x)=T{x)^,f{x), R{x,y)^A-\x,y)-5{x-y)A^\ 
and A — TiXy/A^/^V {V is the system volume, and we 
put the lattice constant a = 1) is the inverse of the DOS 
(mean level spacing) at e = 0. F^\j] is independent of 
the massive modes, F^^ is the purely massive mode, and 
^ is the mixing term. Using the cumulant expansion 
and integrations of the massive modes we obtain Fi ^ 



f['\j] + {fI 



(I)\ 



where 



(i)\ 



R{x,y) 



xy 



X [str T{y)T{x)sti T{x)T{y) 

-str T(y)T(x)E,str T(a:)r(j/)E,] . (14) 



This calculation can be systematically done by using con- 
traction rules derived in Appendix El We neglected con- 
tributions that can be considered higher order ones. The 
first term in the above equation is also neglected since 
the expansion does not include second order in P [see 
Eq.lUni)]. We obtain 



Fi = ^ f R{x,y)siYQ{x)Q{y) 

^ Jxy 



4 



R{x,y)stif{y)T{x)i:,sivf{x)T{y)^, 



ixy 

'2AV 
'2AF 



strE2g(x) 

sir J{x)kT.^Q{x). 



(15) 



The second term has a double-supertrace form and is 
not present in nonchiral systems. The crucial point is 
that the massive modes were parametrized as in Eq. (|ll|l . 
They have the structure E2 in chiral space. 5Q in a form 
5Q — diag((5gi, (5g2) would give the first term of {f[^^) 
only, which is the case for nonchiral systems. 

Using the gradient expansion, we obtain the final form 
of the cr model 



iire 



2AV 



str (VQ)2 - 



strQS^ 



ttDi 
16AV 



(str QVQS,)' 



(16) 



where we neglected the source term, Q{x) — T{x)Y,zT{x) 
is a 4 X 4 supermatrix, and 



ttD _ X-r^a^(r) 



AV 



(17) 



Due to the relation D = Di J^a^(r) ^ Dir^, the con- 
stant Di is smaller than D by the factor 1/rg and the 
second term in Ea. (|16|l can be neglected. However, it 
can be important when the quantum effect is taken into 
account by the renormalization group method. It is dis- 
cussed in Sec llVI 

The generating function Zi is used only for a single 
Green function. It is straightforward to extend the cal- 
culation to the case of products of Green functions. The 
generating function for the product of the retarded Green 
function (trG(^)(ei)trG(^)(e2)) is defined as 



Z2[J] = / X»(V',V')exp 



(18) 



where Vi ^ &re eight-component supervectors. e = 
diag(ei,e2) is the matrix in "two-point" space. In 
chiral symmetric systems, the identity trG*^^^(e) — 
—tiG^^^—e) holds and the generating function for the 
advanced Green function can be found from Z2. Repeat- 
ing the calculation in a similar way, we find the cr model 
{Z2) = f VQexp{-F2) with 



F2 



nD 
4AV 
ttDi 



32AV 

in 
'2AV 



str (VQ)2 

(strQVQS,)' 
str eEzQ, 



(strAQVQS:^)' 



(19) 



where Q — TTi^T is an 8 x 8 supermatrix and A 
diag(l, —1) in two-point space. 



B. Comparison with other models 

Our derived nonlinear a model is equivalent to the 
models in Refs.0, |23| except for the presence of the 
double trace term. The reason why that term was ab- 
sent in Refs.[T3.l20l| is that the massive mode integration 
was not taken carefully. 

In order to compare our result with the models in 
Refs.j3,E| we use the Q-matrix parametrization 



Q = E4(l-p2)i/2+,p]2^ p 



t 
t 



(20) 



where t is a 2x2 supermatrix. The random flux model 
in Ref. 18] is mapped onto the effective action 



Srf — 



b 
2iuj 



J strVT^^VT- (strT-^Vr)^ 



str T + T- 



(21) 



where T G GL(ri|n). This model is reduced to our model 
by using the parametrization 



T= [i+(l+t2)l/2]2^ 



(22) 



and putting n = 1. The "flavor" degrees of freedom n 
represents different species of electrons and are not im- 
portant for the present problem. We note that different 
notation and conventions are used in this expression. In 
contrast with our definition of supermathematics , the 
definition in Ref. 0] was used in Eq. H21() , which explains 
the difference in appearance between Eas. H16|l and 1)21(1 . 

It is worthwhile to mention the relation of the cou- 
pling constants b and c. The authors in Ref. found 
the relation b ^ c/N where N are the "color" degrees of 
freedom. TV must be large in order to justify the saddle 
point approximation. Thus the second term in Ea. ((21|l . 
is small compared with the first term. This is precisely 
what we found and the correlation length rg corresponds 
to TV. We also note that we neglected the topological 
term coming from the boundary condition [l^ . Such a 
term is expected to be derived in our model by consider- 
ing a finite topological number j/ and it will be discussed 
elsewhere. 

In a similar way, our result is compared with the Gade's 
replica a model based on the sublattice models 



Sg = ^ J ti-\/ {Z + W)S/{Z ~W) 
[tr {WV Z ~ Z\/W)f - 



tiW, 



(23) 



where Z is a matrix with some symmetry and W = {I 
^2 -J 1/2^ The parametrization 



2t(l + t 



2x1/2 



W ^l + 2t^ 



(24) 
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is used to find a formal agreement with our model. We 
note that the Gade's model was obtained by using the 
replica method and the structure of the matrix t is dif- 
ferent from ours. However we show in the following that, 
at least in the perturbative regime, both calculations give 
the same result. It is known that the replica and super- 
symmetry methods give the same perturbative result for 
the same symmetry class. 

The relation of the coupling constants b and c was not 
discussed in R.ef . [T^ . It is not clear what is the large 
parameter in the model to justify the saddle-point ap- 
proximation. It is expected that introduction of such a 
parameter leads to a similar relation just as in other cal- 
culations. 

Both works 0, 0| did not use the Q-matrix repre- 
sentation. It has been used in traditional a models and 
is useful for comparison of models and for formulation 
of perturbative and nonperturbative calculations as we 
demonstrate below. It is also important to find gauge 
invariance of the model. 



III. LEVEL CORRELATION FUNCTIONS 

In this section, we calculate the DOS and TLCF by 
using the nonlinear a models derived in the previous sec- 
tion. We neglect the double-trace term contribution and 
put Di — 0. This is because Di is smaller than D by 
the factor l/rg at the classical level. The effect of the 
double-trace term is discussed in Sec lIVI 

We write down the DOS and TLCF in a functional 
integral form. The DOS is given by 



(Pie)) = 



1 



AAV 



Re / VQ 



str k'SzQix) 



'Ft 



, (25) 



where Fi is given by Ea. (|16|l (we put Di = 0) and Q is 
a 4x4 supermatrix. The TLCF is 

(p(ei)p(e2)) = l[Wiei,e2) + W{ei,~e2) 

+Wi-ei,e2) + Wi-ei,-e2)l (26) 
1 



VQ 



16A2y2 

str fcA2S.Q(y) 



str kAiT.;,Q{x 

F2 



(27) 



where F2 is given by Ea. (|19|l . Q is an 8x8 supermatrix, 
and Ai_2 = (1 ± A)/2. In the following we use the con- 
nected part of the TLCF ((p(ei)p(e2))) = {p{ei)p{€2)) - 
{p{ei)){p{e2)). 



A. Summary of the result 

Before entering into the detailed analysis, we give an 
outline of the derivation and the result for reference. For 



perturbation theory, the Q-matrix is expanded in powers 
of the P matrix: 



Q{x) 



l + iP 
1-iP 



(l + 2iP-2P^ 



•)• (28) 



Correspondingly, the result is expressed by using the ex- 
pansion of the diffusion propagator 



n(<?,e) 



A 



2tt Dq^ 



(29) 



The expansion parameter is 1/g where g — ttEc/A = 
ttD/AL'^ is the dimensionless conductance. It does not 
appear in the zero- mode sector of the propagator {q = 0) 
and the expansion is not justified. Actually, treating the 
zero mode exactly (nonperturbatively) , and neglecting 
other nonzero modes, we can obtain the ergodic result. 
The Q matrix for the zero mode is written as 



(30) 



where T is independent of the spatial coordinate and its 
explicit parametrization is given in the following. In or- 
der to incorporate the zero and nonzero modes into the 
analysis we should use the parametrization 



Q{x) = TQ{x)T. 



(31) 



Q parametrizes the nonzero modes and is expanded in 
powers of the P-matrix as in Ea. (|28|) . The zero mode 
Q — TTizT is treated nonperturbatively so that the er- 
godic result is obtained. This parametrization is reminis- 
cent of the renormalization group calculation (see, e.g., 
Ref.Q) and was used by KM. They considered integra- 
tions of the nonzero modes first and found corrections to 
the ergodic result. For a technical reason the result was 
applicable only to the domain z g where z = 7re/A is 
the scaled energy variable. Here we consider the zero- 
mode integration first and then integrate the nonzero 
modes. This method allows us to consider the domain 
2: 3> 1 discussed by AA. For comparison, we present the 
KM method in Sec lmOl 

The zero-mode model is equivalent to chiral RMT. 
This ergodic limit can be obtained by putting g = 00 in 
the above functional integral form. The result is scaled 
by the mean level spacing A to give 



P2{Z1,Z2) 

where 

pf^(^) = 

KiZi,Z2) = 



A {pie = Az/tt)) 



(32) 



= A2((p(ei = Azi/7r)p(e2 = Az2/7r))) 

= -K^{Z,,Z2), 



(33) 



'-[J^(z) + JUz)], 



^^j—^-[ziJiizi)jQiz2) - Z2Joizi)Ji{z2)]. 

Zf - zi 

(34) 
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This result does not depend on any parameter and is 
universal. It was obtained in Ref.|23 by using the or- 
thogonal polynomial method and in Ref.|23 using the 
supersymmetry method. 

How is it changed if we include the nonzero modes? 
If we treat all the modes perturbatively, the result is 
expressed by the diffusion propagator. The expansion 
(|^ is used to give 



((P(£l)p(e2))> 



1 

A 



iRe(^n(<j,, 



, (35) 



^Re^[n2(g,(ei+e2)/2) 



+tf(g,(6i-e2)/2)]. 



(36) 



This expression includes the zero-mode and is justified 
for g ^ 1 and z 3> 1. The zero mode contribution gives 
the asymptotic form of the ergodic result as was shown 
in Ref.jil. 

Before discussing the exact treatment of the zero mode 
we must mention the effect of the renormalization of the 
mean level spacing. The quantity A was introduced as 
the mean level spacing at g = oo and z = oo. For tra- 
ditional symmetry classes, it remains unchanged even if 
we include the nonzero modes (finite-i? effect), which is a 
consequence of the particle conservation law. However, 
this is not the case in chiral systems. For finite-g, the 
nonzero modes contribute to A, which means that the 
DOS is renormalized as was discussed in Ref . [T^ . Refer- 
ring to Ea. (|35|l . we define the renormalized mean level 
spacing 



1 

A 



1 

A 



1 



1, 



-Re 



(37) 



Note that the zero mode is excluded in this expression. 
Contributions of the zero mode are totally different from 
those of other modes. The nonzero modes determine the 
macroscopic behavior of the DOS 1/A, while the zero 
mode determines the universal microscopic behavior after 
scaling in terms of A. 

A naive calculation shows that A is divergent in some 
cases and should be renormalized to a finite value using a 
regularization. We are interested in the microscopic be- 
havior after the mean level spacing is scaled out. The 
effect of nonzero modes in the microscopic domain is 
present even after the scaling and we discuss it in the 
following. 



We turn to the main results in this section. We use 
the parametrization H31I) to treat the zero and nonzero 
modes separately. The zero mode is parametrized so that 
the ergodic results (|32|l and H33|l are reproduced and the 
nonzero modes are treated perturbatively. The domain 
z <C 5 was first considered by KM for nonchiral systems 
and we call it KM's domain. Up to second order in 1/g, 
the DOS in the KM's domain is given by 



A{pie = Az/tt)) 

^(2z- + z' — 
8g2 {'''dz^ dz^ 



1 



ad is the momentum integration 

— y (- 



ad 



pf\z). (38) 



(39) 



We used the periodic boundary condition. The TLCF is 



P2{zi,Z2) = A^((p(ei 



Azi/7r)p(e2 = AzaA))) 
d 



ad ( d 



Z2 



ad 



zi 



d 
dzi 



zi 



d 
dz2 



dz2 



K{zi,Z2) 



(40) 



The result was scaled by the renormalized mean level 
spacing H37() . The calculation of the DOS for chiral sys- 
tems has been done in Ref. |14| but the renormalized mean 
level spacing was not introduced. It leads to a different 
conclusion on level statistics as we discuss in Secs lIII CI 
and El 

We now consider the A A domain z 1, g 1. The 
scaled DOS is given by 



Pi{z) ^ 1 



cos2z„, , 1 



2z ' ■ 82; 
where 2?(z) is the spectral determinant 



(41) 



T^{z) = n 



n 



g^(4^2n^)2 



(42) 



The TLCF is 



P2{zuz2) - ^Re E (H^ + nl) + ^I?ilm J] (H+ + n_) + ^I^^Im J] (H^ - H_) 
-t--^ \ViV2{VlVZ^ - l)cos2(zi -I-Z2) -l-2?iI?2(2?-2?;^ - l)cos2(zi - Z2)] 

8ziZ2 
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1 



[1 + PiX'al'+I?!^ cos 2(zi + za)] 



1 



2(Z1 + Z2)2 

'— 2 (^1 sin 2zi — 7)2 sin 2z2) 



2(zi - Z2)' 



[1 - ViV2V'iV^^ cos 2(zi - Z2)] 



(43) 



where Via = 2?(2i,2), = T>{{zi ± Z2)/2), and ni.2 = 
n(g,ei,2), n± = n(g, (ei ± e2)/2). The result is ex- 
pressed using the spectral determinant as the AA result 
[3- Equation H41|l can be interpreted as follows. Con- 
sider the asymptotic form of the ergodic result 132|l 



Pi{z) ^ 1 



cos 2z 
2z 



1 

8^2 



(44) 



Then, including the spectral determinant in the oscillat- 
ing term, one finds Ea. (|41|l . Ea. l43|l is more complicated, 
but we can see that the ergodic limit gives the asymptotic 
form of the exact result H33|) . While standard perturba- 
tion theory gives nonoscillating terms, expansions around 
two saddle points ^ are required to get oscillating terms. 

We emphasize that Eqs.(|2HI), EOI), ED, and (O are 
the main results in this section. They have the following 
properties. 

Common domain 1 <C z ^ 5. The KM and AA results 
have a common domain 1 <C z <C g where the asymptotic 
expansion of the Bessel function and the expansion of the 
spectral determinant in z/g can be used. In this domain, 
the DOS and TLCF are approximated as 



P2 (21,^2) 



1 - 



cos 2z 
2z 



1 fld 

^ + i^^cos2z, (45) 



sin(zi - Z2) cos(zi + Z2) 



Zl - Z2 



Zl + Z2 



"^(^1 + ^2) cos(zi + Z2) 



ad , ^ . / 

— (zi - Z2)sm(zi 



Z2 



(46) 



Small z. At small energies, the expansion of the Bessel 
function in z is used in Eqs.JSHl and H4(J|) to give 



. ^ nz Od 



^2(2^1, 22) 



Z1Z2 



1 



ad 
2g2 



(47) 
(48) 



These results show that level repulsion at the origin weak- 
ens, which is consistent with the intuitive picture. 

Unitary limit. Taking z,zi -t- Z2 — > 00, we obtain the 
unitary limit as pi(z) — > 1 and 



i?(zi,Z2) = 1 + 



^2(^1,^2) 
Pl{zi)pi{z2) 

1 + Irc y 
2 ^ 

9V0 



cos2(zi — Z2) 

' 2(Z1-Z2)2 



V 



2(zi 

Zl - 



Z2 



(49) 



This result is consistent with the AA's result j^] for 
the unitary class. We note the relations Il{q, e/2; g) = 
2U{q,e;2g) and V{z/2;g) = V{z]2g). The coefficient 2 
in front of g originates from chiral symmetry. Compared 
our cr model (|19|l with the model for unitary symmetry 
0, we see the size of the Q matrix is doubled due to 
chiral symmetry. 

Zl = Z2. The relation p2{z,z) = —p\[z) holds for 
arbitrary g. It can be used to derive the DOS from the 
TLCF. 



B. Density of states 

1. Perturbative calculation 

Now we go into details of the calculation of the DOS 
(|25|1 . The perturbative calculation for nonzero modes is 
considered using the expansion of the Q-matrix in P as 
Ea. (|28|l . The P matrix is parametrized for the chiral 
unitary class as 



t 
t 



p ib 



(50) 



where a, b are real variables, and cr, p Grassmann ones. 
The measure of this parametrization is normalized to 
unity. We define the average 



i>0 (■■■)<=-''' , 

,r^§ls.iVPf-^ls.P^ (51, 

where F^^^ is the second order part of Fi. Performing 
expansions in P as 



(Pie)) 



-Re 1 



1 

2y 



(strfcP^) 



+ 



2V 



(strfep-*) 



(52) 



and using the contraction rules derived in Appendix 
as Ea. l|A6|l . we obtain the result (jSHJ. 

As we emphasized in the previous subsection, this per- 
turbative calculation of the nonzero modes suggests that 
the mean level spacing A is renormalized as Ea. (|37|l . The 
exact definition of A can be written as 



— / strfcS^g(a;) 



p-Fi[Q] 



(53) 
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Thus effect of the self-interacting diffusion bubble is 
renormalized to the mean level spacing. It corresponds 
to imposing the constraint {Q)fi = S^. In Sec lIII Cl we 
give a detailed analysis using the KM method. 



2. Ergodic limit 

At the ergodic limit g ^ oo, spatial dependence of the 
Q matrix is neglected and the DOS is reduced to the form 

pi{z) = jf PQstr/cE^Qexp f-^strE^Qj . (54) 



Following Rcf . [23 , we parametrize the Q matrix as 
Q = TE^T, T = UToU, 

To = 



Of 
i0B 



U = 



u 
u 



u — exp 



e 

ri 



(55) 



where — tt < Of ^ i^ and < 0s < oo. The measure is 
given by 

VQ = deBdOFd^drj^ 

ZTT 

coshes COS — 1 — isinh^B sin ^i? 
^ (cosh 61b - cos 61^)2 ■ ^ ' 

We note that the compact (noncompact) variable 6f (Ob) 
is used for the fermion-fermion (boson-boson) block [g. 
Substituting this parametrization into Ea. (|54|l and inte- 
grating the Grassmann variables, we find 



pi(z) = 1 + Im / dt dOs I d0F- 
Jz Jo Jo ^ 

x{cosh9B cosOf - i)e^*+(™'^heB-cosei.) 

/OO 
dt{jl{t)~Jl{t)). (57) 

Here we introduced the auxiliary variable t"*" =t-\-iQ and 
assumed z > (i > 0). For the Bessel function, we used 
integral representations 



— 2i 

Jo(z)=Re / d^B^ 



cosh Bb 



for the noncompact variable and 

1 r 

Jq{z) = - / dOpe' 
Jo 



(58) 



(59) 



for the compact variable. J\ is given by J\(z) = —Jq{z). 
Integrating the variable t, we obtain Ea. H32l) . 

The asymptotic form at z 3> 1 is given by Eg . (|44|l . 
This cannot be obtained by standard perturbation the- 
ory which gives only nonoscillating terms, the first and 



third terms in Ea . (|44|l . The oscillating second term can 
be obtained by taking into account two saddle points for 
integrals of 9b.f in Ea. H57|l . In addition to the "stan- 
dard saddle point" {0b, 0f) — (0,0) we have another 
"supersymmetry-breaking saddle point" (0,7r). We note 
that the point (0, 0) corresponds to Q = E^ and (0, tt) to 
Q = —kYiz- This is precisely the idea of the calculation 
in Ref. . Taking into account fluctuations around these 
points, we can obtain the desired result. 

In fact this idea is used to find the correct asymptotics 
of the Bessel function. The noncompact representation 
(|58|l is used for 9b and has the saddle point 6'b = 0. 
The compact representation (I59|l for Of has the saddle 
points Of = 0,7r. Expanding around these saddle points 
respectively, we have 



Jo(z) 



1 - 
1 

8^ 



8z 



cos z 



(60) 



It is interesting to note that the expansion around a single 
saddle point is required for the noncompact representa- 
tion (|58|l and two points for the compact representation 
(jSnil- When Ea. (|^ is deformed to Ea. f^ the single 
point = splits into the two points Of = 0,7r. It can 
be shown by considering the deformation of the integral 
contour used in Ref. [23]. We find 



-2i 



dOe 



iz cosn 



2 r'^ 



2i 



dOc' 



dOe 



-z sinh f 



(61) 



This representation is known as the Hankel function 
^0 = Jo + iNo- Taking the real part, we obtain 



1 r/^ 

Jo{z) ^ - I d6'e'"™"' 



1 



7r/2 



dOe- 



(62) 



This expression is reduced to Ea. (|59|l by changing the 
variable 6 n—Oui the second term. Thus the point in 
the second term is changed to tt. Note that the real part 
of the integral is taken in the noncompact representation 
(|58|l . which gives the second term. 

This method, taking into account a set of nontrivial 
saddle points, is the main idea of the nonperturbative cal- 
culation. It produces the exact result for the unitary class 
and the asymptotic ones for the orthogonal and symplec- 
tic classes 9J. It has been used even for the replica |2^ 
and Keldysh a models. For chiral symmetric sys- 
tems at the ergodic limit, a similar technique has been 
used in Ref. 28] to find the asymptotic result (|44|l . In the 
following, we examine how the effect of nonzero modes is 
incorporated into the asymptotic form. 
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3. Integration of the zero mode 

We write the Q matrix as Eq.JJlJ and use the 
parametrizations 1)28(1 and H50|l for Q, and 1(55(1 for T. 
It is shghtly modified as 

Q{x) = UToUQ{x)UfoU ^ UTqQ{x)%U. (63) 

As a result, Fi becomes independent of the Grassmann 
variables of the zero mode. The preexponential term in 
Eg. ((25(1 is written explicitly using the Grassmann vari- 

I 



ables as 

siv kYjzQ{x) s.iT:kYjzT{)Q(x)TQ 

+2e?7strE,Tog(x)To. (64) 

We neglected contributions that vanish after integrations 
over ^ and "q. The first term does not include those vari- 
ables and we can put — 1 for the integrations. The 
second term is also easily integrated out and we thus have 
(P(e)) = (p(e))i + (p(e))2 where 



(P(^)>i - 

(P(^)>2 = 

/(e,0B,0F) = 



1 



Re / VQJ[Q] 



4AF 
1 

27rA 

-z^ / VQJ[Q]e,^j> 
oz 



str kTizQ 



Re / cLOb 





(cosh cos 6'_F — I + z sinh sin 0i?) 
'^"f ^ n ^^. n Ob, Of), 



ttD 

'iAV 



(cosh 6*3 — cosOp)^ 
Str (VQ)^ 



2Ay 



strS^ToQTo 



(65) 



{p{e-))i gives the perturbative result 1(35(1 without the zero 
mode contribution and is equal to the inverse of the 
renormalized mean level spacing 1/A. (p(e)), includes 
the ergodic result and is nonperturbative. 

We note that the Jacobian J^[Q] contribution exists 
in the present parametrization 1(31(1. It depends on the 
nonzero modes only and can be written as 



J[Q] = exp 



(66) 



This contribution changes the renormalized mean level 
spacing slightly and the scaled DOS pi{z) is not changed 
in our approximation. For this reason, we neglect this 
contribution in the present section. It is treated in Sec II VI 
when we discuss the DOS renormalization. 

Let us turn to the calculation of (p(e))2. The kinetic 
term in Fi does not include the zero mode and is ex- 
panded in powers of P. The second term in Fi (and the 
preexponential term) is expanded as 

strE.ToQTo 
= V {cosh B — cos f) 

+ j [str (A:i^P^) cos 6'f + str {kBP'^) cosh 6*5] 

- j [str (kp^xP^) smOp - istr (kB^xP^) sinh^s] 

+ •••, (67) 

where kp^B = (1 i fc)/2. In the following calculation we 
neglect odd terms in the P-matrix. Their contributions 
give 1/f/^-corrections at most. Another reason to neglect 
them is that they involve a factor sin 6 which goes to zero 



at the saddle points Op — 0,tt and 9b = 0. Using this 
approximation, we find the simplified expression 

I(e,OB,Op) ^^Q^z(\B-XF) U^z^BAB+^zXI.Ap\ ^ 

oz ^ ' 

ApMQ] = j StTkp^B^z[Q-^z], 

(•■■>M„ = jvQ{---)e-^--, 
ttD 



^kin — 



4Ay 



str (VQ)^ 



(68) 



where z — 7re/A, — coshOB, and Xp — cosOp. 
Ap^b[Q] include even powers in P. Introducing the aux- 
iliary variable t, we obtain 

-1 /"OC poo pTT 

(P(e)>2 = —rim / dt / Mb / dOpiXB^F - 1) 
ttA Jz Ja Ja 



xe 



t+(\B-XF) 



l-(i-z) — 
OZ 



AzXb Ab +iz\F Af^ 



' kin 



(69) 



Now the problem is how integrations of the variables 
9b,f are performed. They can be done by noting that 
the variable t in the exponential is shifted to t+ -I- zAb or 
zAp compared with the ergodic limit. For the fermion 
part Op, there is no convergence problem and the Bessel 
function is derived. It is also the case for the boson part 
Ob since the convergence problem does not arise for the 
part including Grassmann variables and the other parts 
are real. The only difference is that we cannot take the 
real part for the expression after integration oi Ob since 
the argument I + zAb includes Grassmann variables. We 
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get the Hankel function Hq = Jq + iNo instead of the 
Bessel function Jq [see Ea. (|61f) ]. However, the imaginary 
part iNq does not contribute to the final result since the 
functional Ab is reduced to a real function in the end. 
This is valid in our approximation keeping contributions 
up to second order in 1/g. Thus we neglect the imaginary 
part and obtain 

x(^Ja{t + zAB)Mt ~ zAf) 



-Jiit + zAB)Jiit- zAf] 



(70) 



The ergodic limit g — oo can be found easily by putting 
^f,b[Q] = 0. We note again that this equation was 
obtained by neglecting contributions including sui9f or 
sinh 9b ■ This approximation is valid up to second order 
in 1/g. It still remains to carry out integrations over the 
nonzero modes. In the following we consider two limiting 
cases. 



4- KM's domain (z <^ g) 

The case z <^ g can be considered using KM's method 
For chiral systems, it was considered in Ref . [T^ . In 
our method, the Bessel functions in Ea. (|70|l are expanded 
in powers of zAf,b ^ 0{z/g) to find 



(P(^)>. 



A 



Re 



1 + ^(^B-A,>,,^Z 



Combining with the perturbative contribution 



(P(e)>i 
we find 

(p(e)) ~ 



1 

A 



Re 



1 



A 



Re 



1+ o 



A 



(71) 



(72) 



+ UiAB^AF)\,^^z 



d 

dz 
d_ 

kin dz ~ dz 



Up to here the DOS is scaled in terms of the bare mean 
level spacing A. We introduce the renormalized mean 
level spacing as 1/A = (p(e))j^. Defining the energy vari- 
able as z = 7re/A, we use the transformation formula for 
a function f{z) 



' ~ dz 



^~-^^^b~Af),,^~z- 



fiS)- (74) 



It is apphed to Eo. lfTS)) to find 

pi{i) = A(p(e = Af/TT)) 
~ Re 



l + lmB-AFn.4,^^ 



(75) 



where 

{{{Ab - AFf))i,in 



{{Ab-Af)X^^~{Ab-Af)1^ 

(76) 



9' 



This is obtained by expanding Q in powers of P and 
using the contraction (|A6p with e = 0. ad is momen- 
tum summation and is given by Eq. H39() . Thus we obtain 
Eq.(|2HI)- 



5. AA 's domain (1 <^ z) 

In the limit 1 ^ z, the asymptotic form of the Bessel 
function H60I) is used to write 



(Pie)) 



1 d , 
9 — —Re— / dt 
2 A dz 



t - z 
4^3 



^iz{AB+AF) 



kin 



.t- Z 



t 



2it+iz{AB-Aip) 



A 



1 



-P(.,l,l)--e-P(.,l,-l 



(77) 



where 



I?(z, As, A_f) 
F{z, \b, Xf) 



izXn 



izXp 



2V 



str kp^ziQ — 



2V 



strfcBE,(Q-S,). (78) 



F(z, 1,1) = Fi does not break supersymmetry, which 
means it does not include the supermatrix k — 
diag(l,— 1). As a result we obtain I?(z,l,l) = 1. On 
the other hand, i^(z,l,— 1) breaks supersymmetry and 
the function I? (z, 1,-1) is not normalized to unity. It 
is calculated as I?(z,l,— 1) ^ 2?(z), where the spectral 
determinant 'D{z) is given by Ea. H42(l . We used the ap- 
proximation of keeping second order in P for F(z, 1,-1). 
We refer to Appendix^for details (see also the following 
paragraph) . 

Ea. H77|l is rewritten in terms of the energy variable 
scaled by the renormalized mean level spacing z = ire/ A. 
We use the formula (|74|l and the difference between A 
and A is expressed by the diffusion propagator Il{q,e). 
It represents the self-interacting diffusion bubble and 
should be canceled out. Actually we have contributions 
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from the function X'(z, 1,— 1) by keeping higher-order 
terms in P. Wc find 



V{z,l,-l) 



VQe 



)(l-f^W(z,l,-l)) 
(79) 



1 + izRc[^n{q,€ 

q^Q 



where F^") denotes nth order part in tlic expansion. The 
second term cancels with a contribution coming from the 
transformation (|74|l . Noting 'D{z,g) = 'D{z,g), where 
g = ttEc/ A and g — ttEc/A, we finaUy obtain the result 
Eq.imi) for 1 < z. 



C. Comparison with the KM's method 

The obtained result (|38|l for the KM's domain differs 
slightly from Eq.(21) in Ref. 14] by the presence of mo- 
mentum integration of the propagator X^^^^o ^(9, 0). As 
we can understand from Eq.(j23), the difference comes 
from the introduction of the renormalized mean level 
spacing [Eq.lESl coincides with Eq.(21) in Ref. [13]. It 
is expressed as a self-interacting diffusion diagrams [it 
can be understood by noting the coordinate representa- 
tion n((j) = n(x, x)] and is renormalized to the mean 
level spacing. In order to make this difference clear, we 
repeat the calculation using the KM's method considered 
in Ref. this method, the nonzero modes are in- 

tegrated out while keeping the zero mode variables. It 
allows us to obtain the renormalized effective zero-mode 
action and is useful to understand how we can introduce 
the renormalized mean level spacing. 

We start from the functional for the DOS with the 
source term 



F 



ttD 
4AV 



str [Vg(a;)]^ 



2AV 



it: J 
'2AV 



str kQ{x)'Ez. 



(80) 



The Q-matrix parametrization 1)31(1 is substituted and 
the nonzero modes Q are expanded in P as Eq.J^HJ- In 
our approximation keeping second order in 1/g, the ex- 
pansion is performed up to fourth order in P. The func- 
tional F consists of four parts: 



F 



Fo + F + Fj + Fj 



(81) 



Fq is the zero mode part Fq = F[Q = TS^T], F the 
nonzero mode part F — F[Q], Fi the mixing part, and 
Fj the source term. They are expanded in P as 



F 
Fi 
Fj 



F 



(0) 



J 



F 



(2) 



J 



(4) 

r 

(3) 



F 



(4) 



(82) 



where F^^"^ denotes the nth order part in P. 

The effective functional is obtained by integrating the 
nonzero modes. We define F^ff as 



(83) 



where the average is performed with respect to F^^^. We 
use the contraction rules derived in Appendix^ Up to 
second order in cumulant expansion. 



lire 
2A 



Fo + Ff + (FD + (Fr ) - ^{{FD) - {{Fr'Fy')) 



^(4)\ 



p(2)2\ 



177,7 

2A 



9#0 



str QT,2 



8A2 



(2)^(2)^ 

^tf(g,e)) (strgS,)' 



9#0 



1 + ^ (5] n(g, e)) ' str fcS.g + ^ (J2 n'('?, e)) str gs.str fcE.g. 



(84) 



g#0 



r 



Since momentum summations potentially involve diver- 
gence, this expansion is somewhat cumbersome. This can 
be clearly seen by considering the KM domain z ^ g. 
Then the energy e in the propagator is neglected in our 
approximation n(g, e) n(g,0) and the effective func- 



tional can be written as 

„(1)2 



Fnff ~ — — 



2A 



1 



8g2 



str gs^ 



ITT J 

" 2A 



(1)2 

1 + 



(strgs,)^ 



8A2 4^2 



str kYizQ 
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+ 



str Q'SzSti kTizQ, 



4A2 4c,2 

where ad is given by Ea. (|39|l and 



,(1) 



E 



(85) 



(86) 



P matrix is given by 



This summation is divergent at d > 2 and we need some 
regularization. Fortunately, and as it should be, the 
quantity a^^"^ can be renormalized to the mean level spac- 
ing by defining the renormalized spacing 



1 

A 



1 

A 



Jl)2 



0(1/5^) 



(87) 



This is nothing but the expression (|37|) at the KM do- 
main, a^^'' corrections come from the average {Q{x)). On 
the other hand the second and fourth term in Ea. H85() 
come from the contraction {{Q{x)Q{y))) and cannot be 
renormalized to A. They give the corrections obtained 
in Eq.JSHl)- 

Thus the KM's method makes the problem of the 
renormalization transparent. The idea of integrating out 
fast variables matches the philosophy of the renormaliza- 
tion. Nevertheless, we did not use this method for the 
reason that it is not convenient for calculations in the AA 
domain z 3> 1. Integrations of zero-mode variables nat- 
urally bring contributions from nontrivial saddle points, 
which is an important idea for nonperturbative calcula- 
tions. 



D. Two- level correlation function 

Now we turn to the calculation of the TLCF (|26|) . Q 
is an 8x8 supermatrix and the explicit parametrization 
is different from the previous case. 

For the standard perturbative calculation, we use the 
expansion Ea. (|28|) . The explicit parametrization of the 



P = 

h = 

tl2 = 



t 
t 



Pi ibi 

c irj 
C id 



t21 



t2l - 



tl2 
<2 



02 cr2 

P2 ih 



c 

irj* 



id* 



(88) 



0-1,2, ^1,2 are real variables, c, d complex variables, and 
the greek symbols denote Grassmann variables. As the 
explicit parametrization implies, ti_2 represent the "chi- 
ral" part and ii2.2i the "unitary" part. Starting from the 
expression (|27|l . we have 

I - ^ j^siTk^lP''{x) + ■■■ 
l-^ j^strkk2P\y)+--- .(89) 

,^1, 22, 2^2) given by Ea. (|Alip is second order in 
P and is the base of the perturbative expansion. The 
contraction rule given by Ea. (|A14|l is used to evaluate 
the above expression. The leading order contribution to 
the connected part comes from the contraction 

(strMiF2(a;)str M2P2(y)) = AU^ix - y, (ei + e2)/2). 

(90) 

Thus we obtain the result (|36|l which is valid at 5 ^ 1 
and Zi^2 = '"■£1,2/^ ^ 1- 

The ergodic limit g ^ 00 was considered in Ref . . 
The Q matrix is parametrized as 



Tch is the chiral part 



(91) 



T, 



(0) 



ch 



COS ■ 



V 



—i sin -J 









cosf 


— z sin ■ 



cos^ 






V 





Uch2 

Uchi 

Mch2 



> 

-I sm ^ 


cos^ y 

Uchl,2 = CXp 



and Tu the unitary part 



9i 
§2 



0-1,2 
Pi, 2 



IF 














idiB 














()2F 














id2B 



(92) 
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y(0) ^ 



Uu = 



/ 


cos ^ 

















cos f 


— ie~"^ sin 


1 







— le"^ sin 


§ cost 




V 


-ie-'^ sin # 














Uui 













Uu2 












V 


Mul 





)• 


Uui = exp ^ 


-r 





Uu2 









ie*'^ sin § 






cos § / 



f7 = 



Uu2 = exp 





i^B 



iyy 
-i77* 



V = 



ipF 

(/^B ; ' 



(93) 



cr, p, ^, and f] are Grassmann variables. The integration ranees of the real variables 9, fi, and are chosen properly 
according to the compact or noncompact parametrization 23] . The measure is given by 



T^n M M A ^ ^ cosh6'iBCOs6'iF - 1 ,„ ao a a ^ cosh 6*23 cos ^a;^ - 1 
DQ ~ dtiiBdaiFda-idpi — - — — --ay2saP2F«o'2«P2- 



47r (cosh^iB — cos6'ii?)^ 



'An (cosh6'2_B — cos6'2f) 



'^^B dipF , ^ , sinh fifi sin 4coshf2B cosI^f 

'X.dilBdilF— — ci?«5 drj drj- 



27r 27r 



(coshi7B — cosVLfY (coshi7B + cosfiF)^ 



(94) 



Using this parametrization, after a laborious calculation, we can obtain Ea. (|33|l (see Ref. 23] for the details). 

The nonperturbative calculation using the parametrization (|31|l can be done in the same way as that of the DOS. 
First we integrate the zero-mode variables. The details are presented in Appendix^ and we find for W 



W = W1 + W2, 









({[ 


TT 


d 


2 dzi 


A 


1 + 


TT 


d 


2dz2 



1 + -{Abi-Afi] 



dti{ti - zi) [Jo(ti + ziABi)Jo{ti - ziAfi) - Mh + ziABi)Ji{ti - ziAfi)]| 



o / dt2it2 - 22) [Joih + Z2AB2)Ja{t2 ^ Z2AF2) - Mh + ^2^52)^1(^2 - 22^F2)] 

2 dz2 J^^ 



W2 

/(z,s) 



I dsids2-;-^^^-\-^I{z,s), 



1 

A2 

2 



■{si 



TT 



— <^ZlZ2(si - S2 + Cl)(si + S2 - Dl)Jo yZlSl + Zl 

C2-D2 



Jo Z1S2 - Zl 



Ci+Di 



X(si - S2 + C2)(si + S2 - D2)Jq Z2S1 + Z2 



Jo Z2S2 - Z2 



C2 + D2 



kin 



kin 



(95) 



where si = coshils, S2 — cos 51 i?, and 
1 

2V 
1 



ABia = / strfcBAi,2S2(Q- S^), 



^Fl,2 — — 



2V 



str kF^i,2^z{Q - ^2 



Ci,2 = -[siAb + s2Af±{Aba + Afa)], 
Di.2 = ^[-siAb + s2Af±{-Aba + Afa)] 



A 



F.B 



= A 



F,B1 



A 



F,B2, 



A 



F,BA 



= A 



F,B1 



A 



F,B2- 



(96) 



We neglected odd terms in P as before. Wi includes 
the perturbative contributions H3t)|) and W2 includes the 
ergodic result (jS^ . 

In the KM domain zi,2 ^ .9, the expansion in zA can 
be used. The integrations of the nonzero modes are cal- 
culated up to second order in 1/g. Introducing the renor- 
malized mean level spacing we have 



(p(ei)p(e2)) ^ 



ad ttd f d d \ 
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+ 



+ 2z z —— + z^—'] 
dzi dzi dz2 ^ dz2 J 



(0)^ ^ (0)/' ^ 
Pi (^i)Pi (^2) 



K\zi,Z2) 



(97) 



where zi^2 = 7rei,2/A, and A is given by Ea. lPTTjl . Sub- 
tracting the disconnected part, we derive Ea. (|40|l . 

The AA's domain 1 <C zi,2, 1 ^ 5 is considered using 
the asymptotic form of the Bessel function (|6()|l . The 
details are presented in Appcndix^l From Wi we obtain 
the first part 



9V0 



sin 2zi 
2zi 

sin 2z2 



Viim J2 (n+ + n_) 



9V0 



2Z2 

C0s2(zi + Z2 
8ziZ2 

cos2(zi — Z2) 

8ziZ2 



V2im J2 (n+ - n^) 

9V0 



■ViV2{V%VZ - 1) 



ViV2iVlV-^ - I] 



(98) 



where Vi^2 = 2?(2i,2), V± = V{{zi ± Z2)/2), and n± = 
n((7, (ei ± £2)72). The first term represents the purely 
perturbative contribution. The second connected part 
W2 is calculated in the same way. We obtain 

1 r 1 - c^^^^^i -^''^ViV2V'^ VY^ 



i (e^^^iPi 



2(zi + Z2)2 
e2*^2l?2) 



(99) 



The derived expressions are written in terms of the un- 
renormalized quantity A and we must carry out rescaling 
in terms of A. Additional contributions coming from the 
rescaling should cancel out with terms we did not show 
explicitly here. This situation is the same as the DOS 
case and we finally arrive at Ea. (|43|l . 



IV. THE DOUBLE-TRACE TERM AND THE 
DOS RENORMALIZATION 



In the previous section, we neglected the second term 
in Eas. (|16(l and (|19|1 . This double-trace term includes 
nonzero modes only and changes the perturbative result. 
It appears only in systems with chiral symmetry and we 
therefore concentrate on the DOS. 

At second order in P, we have instead of Eq. H51|l 



lire 
'AV 



str 



(100) 



The presence of the second term modifies the contrac- 
tion rule as Ea. HA9|l . In this case perturbation theory is 
formulated by expansions in 11, Ea. H29() . and 



n2(g,e) 



A 



(101) 



The corresponding expansion parameters are 1/5 oc 1/D 
and gi/g^ cx Di/D^. 

The perturbative expansion gives the DOS 



(Pie)) 



A 



Re 



i + ^n2(<z,6) + l(^n(g,6 



(102) 



The new propagator 112 contributes to the DOS at one- 
loop order. The renormalized mean level spacing is de- 
fined as the inverse of Eq. (|102|) excluding zero-mode con- 
tributions. Actually this result was derived by Gade us- 
ing the renormalization group method jl7l | . In our model 
(|16|l . following the calculation in Ref.|g], we can obtain 
the same renormalization group equations at one-loop 
order as 



db 



eb, (3, 



dc 



c 



-,(103) 



where b ^ 1/g and c ~ l/si- We used e expansion, 
e — d ^ 2, and /i is the renormalization scale. Pb,c are 
the beta functions for b and c. ^ is the zeta function for 
the wavefunction renormalization and corresponds to the 
result in Ea. H102|l . These equations imply a divergence 
of the DOS and delocalization of eigenstates in two di- 
mension. Thus the presence of the double-trace term 
changes the behavior of the DOS significantly. We note 
that the renormalization procedure produces the double- 
trace term even if we start the analysis from a model 
without that term. The quantum effect in two dimen- 
sion increases the coupling constant c. 

We consider the scaled DOS pi{z) to examine how the 
double-trace term contributes to the result. The per- 
turbative result (|102|l is renormalized to the mean level 
spacing. In a similar way as the calculation in the previ- 
ous section wc find in the KM domain z <C g 



Piiz) 



ad 
8g2 



2z— + z' 
dz 



Od 
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p'^\z). 



(104) 
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In the A A domain 1 <§i z, the spectral determinant is 
modified as Ea. (|A10|l . Subtracting the renormahzation 
effect, we obtain 



V{z) 



n 



g>0,g2^0 



(Dq 



,2^2 



(105) 



which is consistent with Ea. p()4(l . 

Finally we mention the Jacobian contribution in 
Ea. (|66|l . It includes a term second order in P and changes 
the contraction rules. Since this term is similar to the last 
term in Ea. HlUU| it can be easily incorporated into the 
contraction rules by the replacement 



n2(g,e)^n2((Z,e)-n2(q,e). 



(106) 



Thus this Jacobian contribution is always subleading 
compared to the propagator 112. We also note that this 
contributes only to A and not to the scaled DOS pi{z) 
in our approximation. 



V. DISCUSSIONS AND CONCLUSIONS 

We have studied disordered systems with chiral uni- 
tary symmetry. Using a chiral symmetric random matrix 
model we derived the nonlinear tr models (|16|) and (|19|1 . 
We demonstrated that they are equivalent to related chi- 
ral symmetric models. Using the a models, we calculated 
the level correlation functions. We exploited the nonper- 
turbative methods developed by Kravtsov and Mirlin, 
and Andreev and Altshuler for the traditional classes. 

The equivalence of models shows the universality of 
disordered systems. Our derived a models are applicable 
to models treated in Refs.p^ flM l20j. The double-trace 
term was not derived in Ref.|2^. This is because the 
massive mode integration was not considered carefully. 

For the calculation of the DOS and TLCF, we stressed 
the need for the renormahzation of the mean level spac- 
ing. This renormahzation is absent in traditional nonchi- 
ral systems. After separating the renormahzation effect, 
we found the results (|38|l and H40I) in the KM domain, 
and H41I) and (|43|1 in the AA domain. It is interesting to 
note that the results in the AA domain are expressed us- 
ing the spectral determinant. It contributes to oscillating 
terms only, in a similar way as for the traditional classes. 
Thus we conclude that the singularity of the form fac- 
tor at the Heisenberg time is washed out due to fmite-g 
effects 0. 

Our formulation of the perturbative and nonpertur- 
bative calculations can be useful not only for the level 
correlation functions but also for the conductance and 
other quantities. In the present work we concentrated 
on the level correlation functions. In Ref.|23|, the same 
quantities were calculated perturbatively. The different 
result obtained there is due to another parametrization 



of the Q matrix. Additional contributions coming from 
the integration measure would give the correct result. In 
Ref . [l3| ■ the DOS in the KM domain was calculated from 
the model derived in Ref. ^20). The result was scaled in 
terms of the bare mean level spacing A, and the renor- 
malized mean level spacing A was not introduced. This 
leads us to a different conclusion on level statistics as we 
mention below. 

Let us discuss the importance of introducing the renor- 
malized mean level spacing. There are numerous works 
on the behavior of the DOS at the origin e = 0. The main 
question is whether the DOS diverges or not, and ana- 
lytically it has been considered using perturbation the- 
ory at weak disorder g ^ 1. On the other hand, chiral 
RMT, which corresponds to the model at g = oo, pre- 
dicts the vanishing of the DOS at the origin z — 0. This 
is not a contradiction and indicates the importance of 
scaling. The macroscopic behavior is determined by the 
nonzero modes and a divergence of 1/A was reported 
in Ref. The zero mode has nothing to do with this 
behavior and determines the universal behavior at the 
microscopic scale. It can be seen by scaling the energy 
variable e in terms of the mean level spacing. 

Generally speaking, the behavior at the macroscopic 
scale depends on the model. From a field theoretical 
point of view, the divergence can be renormalized to 
the mean level spacing and a definite conclusion as to 
whether it is a real divergence or not can be obtained 
by referring to other approaches such as numerical sim- 
ulations. Our result relies on perturbation theory and 
the divergence may be cut off somewhere before the ori- 
gin. This crossover to the universal microscopic domain 
is highly nonperturbative. Since a high resolution is re- 
quired, it may be hard to see such a crossover numerically. 

From the viewpoint of level statistics, the DOS must 
be scaled (renormalized) to unity at all energies to find 
the universality. This unfolding procedure cannot be ap- 
plied to the present chiral case because the DOS itself 
has the universal fine structure (oscillations due to level 
repulsion) at the origin. For this reason, we use A, the 
(inverse) DOS at z = cxd (e = 0), for scaling in the er- 
godic regime. It is modified by finite-g effects and we use 
A to see the microscopic domain closely. 

Thus using the renormalized mean level spacing, we 
can separate problems at both scales. The effects of 
nonzero modes (finite-g effect) cannot be scaled out com- 
pletely in the microscopic domain and deviations from 
the universal behavior are obtained as we have shown 
in the present work. Such an example can be found in 
Ref. 16]. The generalized random matrix model was used 
there and it was found that the quantity A is different 
from our result. However, after scaling in terms of the 
nonuniversal quantity A, we can find the complete agree- 
ment up to (mite-g corrections. This demonstration of 
"universal deviation" justifies the introduction of A. 

The double trace term contribution is small at the clas- 
sical level because the coupling constant is small com- 
pared with that in the diffusion propagator. However, 
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quantum effects affect tliis coupling and the contribution 
becomes important in some cases. It significantly affects 
the DOS renormalization and a diverging DOS was found 
m Concerning level statistics, this term modifies 

the spectral determinant as Ea. ljA10|) . 

Our calculation is only for the chiral unitary class. The 
other chiral classes, chiral orthogonal and symplectic, can 
be calculated in the same way. The problem is that the 
proper parametrization of the zero mode has not been 
found. However the KM's domain can be considered 
without knowing the zero mode parametrization as was 
done in Rcf. 14]. The obtained result is valid only at first 
order in \/g. Repeating the same calculation up to the 
next order and introducing the renormalized mean level 
spacing, we found the same form as Eq.JSHI)- The coeffi- 
cient of the second term in Eq. (|38|l is changed but with 
the same sign for all the classes. This result also holds 
for the TLCF (|4()(l . We thus obtain the same conclusion 
as KM, namely, the weakening of level repulsion [This 
can be seen, e.g., in Eq.(0HJ)]. The authors in Ref.Q 
drew a different conclusion by looking at the first order 
correction to the mean level spacing. It is renormalized 
to the mean level spacing and should be applied to the 
DOS behavior and not to level repulsion. 

As an interesting application, we mention a related 
work in Rcf. 29] . For traditional nonchiral systems, the 
authors in Rcf. i;|j pointed out that the AA's result is 
related to the Calogero-Sutherland model at finite tem- 
perature. It is shown in Ref.'s^ that this model is equal 
to the generalized random matrix model proposed in 
Ref . [31I] . In this problem the nonlinear a model is mod- 
ified due to power-law correlations of random matrices 
ll.'ij ] and the diffusion propagator and spectral determi- 
nant are modified. As a result agreement with the result 
in Ref.jSlj] was found and a conjecture to more general 
cases was made. We expect this holds also for chiral sys- 
tems and the result is presented in Ref.TS]. 

Another future problem is the wavefunction statistics. 
For traditional classes, it was considered in Ref.js^ using 
the KM method. It will be interesting to see how this 
result is modified in the chiral symmetric case. 



~Xb I strfcs^'^ 



(Al) 



where kF,B = (1 ± k)/2 and z = ire/ A. The P ma- 
trix is a 4x4 supermatrix including nonzero modes and 
is given by Ea. (|50|l . Since this functional breaks super- 
symmetry for Xp As, the function X'i(z, As, Aj?) — 

J VQ exp{— F^^^ ) is not normalized to unity. We calcu- 
late this function and derive the contraction rules. 



Using the explicit parametrization l)50|l. we write 



(2) 



as 



P 

a 
\b 



(9) 



G = 



diag(n((j, eA+), n((j, eX+),U{q, eAi^), n(q, eAs)), 

(A2) 



where A+ = (As -I- Ai?)/2 and the diffusion propagator is 
given by Ea. H29(l . Then the functional integral is given 
by 



■Di{z,Xb,Xf) - j I?(V;,V)e-^i 
"n(g,eAs)n(q,eAi.)T'/' 



n 



n 



n2(g,eA+) 

{Dq^ - ieX+f 
{Dq^ ~ ieXB){Dq^ - ieXpY 



(A3) 
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APPENDIX A: CONTRACTION RULES 



Calculation for Fi 



Consider the functional 



F^'' = ^/str(VF)^-^A. /strfc.P^ 



Since V'('?) ^-nd 'ip{—q) are not independent of each other, 
the square root appears in T>. We used the periodic 
boundary condition and qi = 2Trni/L, rii is integer. 

Using the result we obtain the contraction rules for the 
matrix P as 
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(str AP{x)BP{y)) 



(str AP(x) strBP(y)) 



(str kiAstr kjB — str fc^ylS^str kjBYiz 



i,j = F,B ^ 



(A5) 



r 



where A and B are arbitrary supermatrices and (• • •) = 
V-\z, \b.Xf) J VQ{- ■ ■) exp(-Ff )). 

We are mainly interested in the case {Xb,Xf) = (1, 1), 
and (1,-1). The first case (1,1) corresponds to standard 
perturbation theory. The free energy does not break su- 
persymmetry and we find I'(z, 1, 1) = 1 and 



{sii AP{x) BP (y)) = 



(strAP(a;)strBP(y)) = 



-U{x-y,e) 

X (str Astr B - str AS^str BY.^ 
+str 74E2,str BT,^ 
-str AY^y str BY, y), 

hj{x-y,e) 

xstr {AB - AY^BYz 
+AY^BY^ - AYyBYy). (A6) 



In the case (1,-1), supersymmetry is broken and this is 
used for calculations in the AA's domain 1 <C 2. The 
function V is given by 



Pi(z,l,-1) 



n 



n 



(AT) 



Effect of the double-trace term 



We consider the effect of the double-trace term. The 
second term of Ea. (|100|) is included in Ea. (|Al|) . In this 
case, the matrix G in Ea. ljA2|l is replaced by 



G{q) 



diag(n+,n+,cnF,cns) 

/ 
000 



A 



1 -i 
V -1 



n 



c 



27r Dq^ 



ieXp.B, 
1 



1 



(A8) 



As a result, V and contraction rules are modified in the 
following way. The contraction for {Xb,Xf) = (1,1) is 
expressed as 



{str AP{x) BP {y)) 



1 



n(x - y, e) (str ^str B ~ str AE^str BY^ + str ^E^str BY^ - str AS^str BYy) 
in2(x - y, e)str AYxBY^, 



(str AP{x)stT BP{y)) 



-n(x - y, e)str (AB - AY.BY^ + AY^BY^ - AYyBYy 
- -112(2; — y, e)str AY^str BYx, 



(A9) 



r 



where the propagator 112 in momentum space is given by 
(Uni}. For (As, Xf) = (1,-1), ISlll is replaced by 



Calculation for F2 



n 



{Dq 



2\2 



Consider 



(AlO) 



Pf )(zi,Z2,Z3,Z4) = ^J str(VP)2 
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J strkpA^P^ -'-y j strksA^P^ (All) 

where Ai^2 = (1 ± A)/2 and ^1^2, 3, 4 = 7rei_2,3,4/A. The 
P matrix is an 8x8 supermatrix and is parametrized as 



Ea. (|88|) . This case is considered in the same way as the 
(2) 

case of ' . We neglect the double-trace term contribu- 
tion for simplicity. The result is expressed for the func- 
tional integral as 



1 2 



V (z z) - n + 

J^ch[Zi,Z2 — II — — ^ ■ — „ : — r, 

,>o ,Vo W-^e^){Dq^-^e2) 

T>u{zi,Z2,Z3,Zi) = [[ — — — — . (A12) 

Dg^ - f (ei -f £3) Dg^ - ^(£2 + 64) 



For example, X'2(-2:i, ^i, Z2, ^2) = T^iizi), 'C>2{zi, zi, -Z2, Z2) = 'Di{z2), and X'2(-2i, zi, -2:2, Z2) = 
I?i(zi)I?i(z2)2'i[(2^i + 2;2)/2]X'f^[(2;i — Z2)/2]. The contraction rule is given by 



{str AP{x)BPiy)) = \ ^ ^ n„,y5 (x - y) 



4 

a,P=F,B i,i=l,2 



X (str fca Ai^str kpAjB — str fc^AiSz^str kpAjY^zB + str fc^ AiEa;ylstr kpAjY^^B — str fc^AiSj^^str kfjAjYjyB) 
{sir AP{x)slTBP{y)) = - ^ ^ n,„,/3(a;-y) 



4 



xstr (fcaAi^fc/jAjB - kaAiYjy_AkpAjY.^B + kaAiY^^AkpAjY^^B - kaAiT,yAkpAjT,yB) , (A13) 

where 1X10^73(2;) = Il{x,{Xia + Xjp)/2), and Aif = ci, Am = £2, A2F = £3, = £4- The case (21,^2,23,24) 
(21, zi, 22, 22) corresponds to standard perturbation and we find 

{sti AP{x) BP (y)) ^ ^ n fx - y, ^^^^ j (strAA.strSAj - str^AiE^strSAjE^ 



-l-str AAiS2,stri?AjI]j; - str AA^Sj^str iJA^Sy), 
(strAP(a;)strBP(y)) = J ^ H ( 2; - y, ] str (AA.BAj - AA.E.SA^S, 



4 ^ V 2 

-AA,J:.^BAjJ:^ ~ AAJ^yBAjY^y). (A14) 
I 



APPENDIX B: CALCULATION OF THE 
TWO-LEVEL CORRELATION FUNCTION 

1. Zero mode integration 

In this section we derive Ea. (|95|l by integrating the 
zero-mode variables of the nonperturbative parametriza- 
tion H31|) . As before, the parametrization is slightly mod- 
ified as 

Q{x) = (7uTchri°)g(x)T^°)Tehf/u, (Bl) 

to eliminate the grassmann variables of unitary part in 
¥2- For the pre-exponential term, dependence of the 



grassmann variables on C/u is explicitly written as 

str fcAiSz(3(x)str fcA2Sz(5(x) 

str kAiY.zTc\,Q{x)fc\,stT kA2^zTchQ{x)f^h 

-4Crw*strAiI],rchT^°)g(x)T'^°)Teh 

xstrA2l],Tehr^°)Q(a;)ri")7;h. (B2) 

The neglected terms do not contribute to integration of 
the Grassmann variables. The first term does not include 
the grassmann variables ^ and 77. We can set T^^ = 1 
and have 
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1 



16A2\/2 



VQ 



str fcAii;^(3(x) 



str kA2^,Qiy) 



-F2[Q] 



(B3) 



where Q = TchQTch- It still includes the zero-mode variables of the chiral part Teh- Since the chiral part parametriza- 
tion is the same as that of the DOS, the calculation can be done as in Sec lIII Bl As a result Wi{ei, £2) in Ea. l|95|l is 
obtained. It includes a perturbative part and connected and disconnected parts. 

Next we consider the second contribution which includes only the connected part. It is obtained by integrations of 
^ and 77 as 



W2{e,,e2) = — J ds,ds2 — — j^r^Iizi.2,^i.2,^B,F), 

I{zi,2,si,2,^b,f) = -7^7^ VQVQ,he-^^^\ 
dzi dz2 J 



(B4) 



where Q{x) = T^y^T^^ Q{x)TlP^Tch_, Qch — TchS^Tch, and zi.2 — 7rei^2/A. We examine str eY,zQ{x) to integrate out 
variables in Qch- The expression is simplified if we apply the saddle point approximation we use in the following. At 
the saddle-point we have sinO = and sin 57 = 0. This approximation leads to the reduction 



stTeTizQ{x) — > eistr 



cos f2 -I- A 



cos 57 — A 



cos 6*1 + (cosh 6*15 — co8 9if)pi(7i 'SzQ{x) 

-e2Str ^ — — cos 6*2 + (cosh 6123 - cos 6'2f)/020-2 S^Q(a;). (B5) 
Again we stress that this approximation is justified at second order in 1/g. Substituting this expression, we have 

I{z,.s,ip) = - /l?Qe-^x»/i(zi,s)/2(z2,s), (B6) 



d_ 
dz 



ch 



iz f cos57±A, , ^ ^ , ^ A/ \ 

1 - — / str (cosh6'iB - cos0iF)Pia^T.;,Q(x} 



X exp 



iz f cos 57 ± A X ^ ^, . 
"W / 2 cos9,^zQ{x) 



(B7) 



where _Fkin is the kinetic part in F2 . The variable ip is not included in the integrand in our approximation. Integrations 
of the remaining zero-mode variables are carried out and we find 

Kb\f — 1 



47r J \b - AiF 



1 + iz [ SiX^B - S2\iF + ^(Ci - Di)\iB + ^(Ci -f Di)\iF 



X exp 



iz ( siAjb - S2\iF + ^(Ci - Di)\B + ^(Ci + Di)\iF 



— (Sl - S2 + Ci){Sl + S2 - Di)Jo I ZSl + Z 1 Jo [ZS2 - z 



(B8) 



where XiF = cos 6iF and XtB = cosh 9iB- This result yields W2 in Ea. l|^ . 



2. AA's domain 



We consider Ea. (|95l) in the AA domain 1 ^ zi,2 using the asymptotic form of the Bessel function (|60|l . For Wi, 



Wi{ei,e2) ~ ^( /i(zi)e^^^(^«i+^-i) -f 5i(zi)e2'^i+^^i(^-i-^-i) 
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kin 



'^^^^'^^'^^'^^Vi(^i)/2(^2) + e2^^^-^(-^-^-^-^^)gi (21)32(^2) 



(B9) 



where 



(BIO) 



and F is the supersymmetry breaking functional 



F(zi,Z2,Z^,Zi) 



ttD 
4Ay 



str (vg)2 



-— / StT{zikF + Z2kB)Al^,iQ-^z) 



2V 



str {zskp + ZikB)A2^z{Q - S.).(B11) 



The condition zi = Z2 and z^ = 24 recovers supersymme- 
try. As before we expand the nonzero modes Q in terms 
of the P matrix and use the contraction rules derived in 
Appendix El The first term in Ea. l)B9|l does not break 
supersymmetry [-^(^1, zi, Z2, z^) = F2] and is nothing but 
the purely perturbative contribution. Its connected part 
gives the first term in Ea. H98|l . The second (third) term 
in Ea. (jB9|) gives the second (third) term in Ea. (|98|l . For 



the leading order contribution, we use 
j r)Qe--^(-^i'^i'^^'^^) - Pi, 

j •Dge--^(^i'^i - V2, 

(strMip2)^,„,(^^^^^_^^^^^) - -25](n+-n_). 

(B12) 



For the last term in Ea. (jB9|) . we have 

j x>Qe-'^(-^i'^i'-^^'^^) - V{DiV\VZ^. (B13) 

The disconnected part is included in this contribution 
and is subtracted to give the fourth and fifth terms in 
Eq.(EHI). 

The purely connected part W2 is calculated using the 
asymptotic form of the Bessel function. We obtain 



W2(ei,e2) 



(*) 



(B14) 



where (*) denotes the complex conjugate of the preceding term. Integrations of si^2 are evaluated to find the asymptotic 
form 



W2(ei,e2) 



1 



4A2 



VQ 



1 



'F(zi,Zi,-Z2,-Z2) _ 2i{zi-Z2)-F{-Zi,Zi,Z2,-Z2) 



+ 



1 



{zi + Z2Y 



{Z, - Z2)2 

F (zi ,Zl ,Z2 ,Z2) _j_ ^2i(zi+Z2)~F{-Zl,Zl,-Z2.Z2) 



+ (^1,2 -Zl,2) 



+- 



'-1 - ^2 



g2izi-F(-zi,zi, 22,22) g2i2i-F(-zi,zi,-Z2,-Z2) 



+e 



— 2iz2-F{zi,Zi,Z2.-Z2) I —2i22—F( — Zi,-Zi, 22,-22) 



Ul,2 -Z1.2, 



-ZI.2) 



(B15) 



Finally, keeping second order in P for the functional F 
and using the formula (|A12p . we obtain the result (|99|l . 
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